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Abstract 

The dynamics of the delay logistic equation with complex parameters and arbitrary 
complex initial conditions is investigated. The analysis of the local stability of this 
difference equation has been carried out. We further exhibit several interesting 
characteristics of the solutions of this equation, using computations, which does 
not arise when we consider the same equation with positive real parameters and 
initial conditions. Some of the interesting observations led us to pose some open 
problems and conjectures regarding chaotic and higher order periodic solutions 
and global asymptotic convergence of the delay logistic equation. It is our hope 
that these observations of this complex difference equation would certainly be an 
interesting addition to the present art of research in rational difference equations 
in understanding the behaviour in the complex domain. 
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1 Introduction and Preliminaries 

Consider the delay logistic difference equation 


= ,n = 0,l,... (1) 

t + PZn-l 

where the parameter a is complex number and the initial conditions Z-i and zq are 
arbitrary complex numbers. 

The same difference equation including some variations of the this are studied when 
the parameter a and initial conditions are non-negative real numbers, Eq.Q was inves¬ 
tigated in P] and [H]. In this present article it is an attempt to understand the same in 
the complex plane. 

The set of initial conditions z_i, Zq eC for which the solution of Eq. Q is well defined 
for all n > 0 is called the good set of initial conditions or the domain of definition. It is 
the compliment of the forbidden set of Eq. Q for which the solution is not well defined 
for some n > 0. It is really hard to find out either the good set or forbidden set for 
the second and higher order rational difference equations due to the lack of an explicit 
form of the solutions. Eor the rest of the sequel it is assumed that the initial conditions 
belong to the good set [Hj . 

Our goal is to investigate the character of the solutions of Eq. 0 when the parameters 
are complex and the initial conditions are arbitrary complex numbers in the domain of 
definition. 

Here, we review some results which will be useful in our investigation of the behavior of 
solutions of the difference equation 0 

Let / : —7- D where D C C be a continuously differentiable function. Then for any 

pair of initial conditions zo,Z-i G D, the difference equation 

Zn+l ~ fiZn,Zn—l) , ( 2 ) 


with initial conditions Z-i, Zq G D. 

Then for any initial value, the difference equation (1) will have a unique solution 
A point 2 G D is called equilibrium point of Eq. Q if 

f{z,z) = z. 

The linearized equation of Eq. Q about the equilibrium ^ is the linear difference equation 

Zn+l = aoZn + aiZn-l , u = 0, 1, . . . (3) 
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where for i = 0 and 1. 


of _x 
a-i = ■^{z,z). 

OUi 

The characteristic equation of Eq. ([^ is the equation 

— oqA — ai = 0. 


( 4 ) 


The following are the briefings of the linearized stability criterions which are useful in 
determining the local stability character of the equilibrium ^ of Eq. @), m 

Let ^ be an equilibrium of the difference equation z^^i = f{zn,Zn_i). 

• The equilibrium ^ of Eq. ( 2 ) is called locally stable if for every e > 0, there 
exists a 5 > 0 such that for every zq and Z-i G C with \zo — z\ + \z-i — A| < (5 we 
have \zn — z\ < e for all n > — 1 . 

• The equilibrium z of Eq. ( 2 ) is called locally stable if it is locally stable and if 
there exist a 7 > 0 such that for every Zq and z_i G C with |zo — z| + |z_i — z| < 7 
we have lim^^oo z„ = z. 

• The equilibrium z of Eq. (2) is called global attractor if for every zq and z_i G 
C, we have lim^^oo 

• The equilibrium of equation Eq. ( 2 ) is called globally asymptotically stable/fit 
is stable and is a global attractor. 

• The equilibrium z of Eq. (2) is called unstable if it is not stable. 

• The equilibrium z of Eq. (2) is called source or repeller if there exists r > 0 
such that for every zq and z_i G C with |zo — z| + |z_i — z| < r we have |z„ — z| > r. 
Clearly a source is an unstable equilibrium. 


2 Local Stability of the Equilibriums and Bounded¬ 
ness 

In this section we establish the local stability character of the equilibria of Eq. Q when 
the parameters a and (5 are considered to be complex numbers with the initial conditions 
Zq and z_i are arbitrary complex numbers. 

The equilibrium points of Eq. Q are the solutions of the equation 

az 

1 + /3z 
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Eq.(j^ has two equilibria points zi ^2 = 0, respectively. The linearized equation of 
the rational difference equationQ with respect to the equilibrium point = 0 is 

^n+l ^ 0, 1, . . . (5) 

with associated characteristic equation 

- aA = 0. (6) 

The following result gives the local asymptotic stability of the equilibrium zi. 
Theorem 2.1. The equilibriums Zi = 0 of Eq. 0 is 

locally asymptotically stable if and only if 

lal < 1 


unstable if and only if 


a| < 1 


and non-hyperbolic if and only if 

\a\ = 1 

Proof. The characteristic equation of the equilibrium as already mentioned above is 
A^ — «A = 0. The zeros of this polynomials are 0 and a. Therefore the result follows 
trivially by Loeal Stability Theorem. 

□ 


Theorem 2.2. The equilibriums Z 2 = o/Eg. is 
locally asymptotically stable if and only if 


unstable if and only if 


1 


Proof. The characteristic equation of the equilibrium ;22 is — A + = 0 . The 

zeros of this polynomials are a+y/{4-3a ^^ equilibrium Z 2 is locally 

asymptotically stable if and only if the modulus of the zeros of the characteristic equation 
are less than 1 . Simple algebraic calculation of these two inequalities reduced to the 
condition | < |q;| < |. In the similar fashion the equilibrium is unstable if and only 
if the modulus of the zeros are greater than 1. Consequently, the condition reduces to 
0 < |q;| < |. Hence the theorem is proved. 

□ 

Now we would like to try to find open ball 5(0, e) e C such that if Zn G 5(0, e) and 
Zn-i e 5(0, e) then Zn+i G 5(0, e) for all n > 0. In other words, if the initial values 
Zq and z_i belong to 5(0, e) then the solution generated by the difference equation Q 
would essentially be within the open ball 5(0, e). 

Theorem 2.3. For every e > 0 and for any eomplex parameter |«| < 1 and |/3| < 1, if 
Zn and Zn-i G 5(0, e) then Zn+i G 5(0, e) provided, 

1 

0 <e< - 



Proof. Let {zn} be a solution of the equation Eq.([^ Let e > 0 be any arbitrary real 
number. Consider Zn-i G 5(0, e). We need to find out an e such that Zn+i G 5(0, e) 
for all n. It is follows from the Eq.([^ that for any e > 0, using Triangle inequality for 
\a\ < 1 and |/5| < 1 


l^n+l I 


CXZn 

1 + P^n-l 


\a\ e 

~ 1 - l/^l e 


In order to ensure that \zn+i\ < e, it is needed to be 


\a\ e 

l-\(d\e 


< 1 la 


1 < — |/d| e e < 


1 — |a 


Therefore the required is followed. 


□ 


3 Periodic Solutions 

The global periodicity and the existence of solutions that converge to periodic solutions 
of the difference equation ([^ are adumbrated in this section. 

A solution {zn}n of a difference equation is said to be globally periodic of period t 
if Zn+t = for any given initial conditions, solution {zn}n is said to be periodie with 
prime period p if p is the smallest positive integer having this property. 
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3.1 Existence and Local Stability of Prime Period Two Solu¬ 
tions 

Here we find out prime period two solutions followed by the local stability analysis of 
them. 


3.1.1 Existence of Prime Period Two Solutions 


Let ..., 0, (p, 'ip,---, (p 7^ ^ be a prime period two solution of the dif¬ 


ference equation Q . Then 


a'lp 


and Ip = 


acf) 


This two equa- 


__ 1+/30 ^ l-\-l3'ip 

tions lead to the set of solutions (prime period two) except the equilibriums 

(-0.5-0.5a);3-0.5-y/(l+(-2-3a)a);92 (-0.5-0.5a)/9+0.5^(l+(-2-3a)a)/92 

^2 ; 'ip ^2 


0 and ^ 


as 




/32 


and 


^ ^ (-0.5-0.5Q;)/3+0.5Y^(l+(-2-3Q)a)/32 ^ (-0.5-0.5Q;)/3-0.5Y^(l+(-2-3Q)a)/3 

0 ^2 ) 'ip fl2 


/32 


3.1.2 Local Stability of Prime Period Two Solutions 

Let - - - ,0, 'ip, 0, 0,.. 0 ^ 0 he 0 . prime period two solution of the Pielou’s equation 

Q. We set 

^n—1 

Then the equivalent form of the delay Logistic difference equation Q is 


^n +1 

_ ChVn 

^n+1 1 I /Q 

1 + pUn 

Let T be the map on (0, oo) x (0, oo) to itself defined by 



Then ( ^ 


is a fixed point of T^, the second iterate of T. 


rp2 


2 U 


rji2 


2 ( U 




av 


\ 



f 9 {u,v) \ 

V h0i,v) ) 


V 
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OLV 


where g{u, v) = and h{u, v) = . Clearly the two cycle is locally asymptotically 


stable when the eigenvalues of the Jacobian matrix Jt^, evaluated at 
the unit disk. We have, 


V’ 


he inside 


1 / 7^2 


4> 

xp 






where (i+W)2 i+^4> 

^(rh 7 /d =_3rirl ^(fh 7 /d =__|_ 

(1+/3(/))2(1+/3'0) (1+/30)(1+/3'0)2 ' (l+/30)(l+/3'0) 

Now, set 

Sg,, + “l+wtO 


(1 + /50)2 


^ = (TT/^oWW 

In particular for the prime period 2 solution. 




(-0.5-0.5a)/9-0.5^(l+(-2-3a)a)/92 




(-0.5-0.5a)/9+0.5^(l+(-2-3a)a)/92 


a 


X = 


(-0.5-0.5a)/3+0.5^(l+(-2-3a)a)/32 ^ “C+ 


(-0.5-0.5a)^-0.5V(l + (-2.-3.a)a)^2 ' 


^ (-0.5-0.5a)^+0.5V(l + (-2-3a)a)^2 y 


and 


(-0.5-0.5a)^-0.5-y/(l+(-2-3a)a)^2 


(^(-0.5 - 0.5«)/5 + 0.5y/(l + (-2 - 3«)«)/52y 


A =-^— ; - 

P2 (-x (-0.5-0.5a)/3-0.5-y/(l+(-2-3a)a)/32\ / (-0.5-0.5a)/3+0.5-y/(l+(-2-3a)a)/32 

p yi--\ ^ j ^ 


and for the prime period 2 solution. 
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(-0.5-0.5a)/9-0.5-v/(l+(-2-3Q)a)/32 “I 1+ 

Ot 1-1-h 


(-0.5-0.5a)^+0.5A/(l + (-2.-3.a)a)^2 


X 


(I (-0.5-0.5a)/3-0.5\/(l + (-2-3a)o;)^2 \ ^ 

I ) 


and 

A = 


^ (-0.5-0.5a);9+0.5-y/(l+(-2-3a)a);92 


((-0.5 - 0.5a)/5 - 0.5^(1 + (-2 - 3a)a)/32y 


A (-0.5-0.5Q)/3+0.5^/(l+(-2.-3.Q)a)/32\ A (-0.5-0.5a)/3-0.5x/(l+(-2-3Q)a)/32 

f* (1 +- -e - ) (1 +-F- 


Then it follows from the Linearized Stability Theorem that both the eigenvalues of the 
Jx 2 ^ lie inside the unit disk if and only if|x|<l + |A|<2. 

In particular, for a = z and (3 = 2 + 3z the prime period 2 solutions is 

{<p —>■ —0.294567 + 0.313317z,V' ^ —0.0900486 — 0.236394z}. The periodic trajectory is 
shown in Fig.l of which the \x\ — 1.38112 and |A| = 0.689678. By the Linear Stability 
theorem (|x| < 1 + |A| < 2) the prime period 2 solution is locally asymptotically stable. 



Figure 1: Periodic Cycle solution of period 3 for a = z and /? = 2 + 3b 

For a = 1 + i and /? = 2 + 3i the prime period 2 solutions is 0 —0.168166 + 0.5344112, 

xp —7- —0.370295 — 0.226718b The periodic trajectory is shown in Fig. 2 of which the 
lx I = 1.5 and |A| = 1.58114. Hence the prime period 2 solution is unstable. 



Figure 2: Periodic Cycle solution of period 3 for a = 2 and /? = 2 + 3b 






















3.2 Higher Order Cycles in Solutions 

Here we shall explore the higher order periodic cycle of the difference equation Q in a 
vivid manner. 

For the parameters a = i and /5 = 2 + 3i of the difference equation, the period 3 cycle 
is the following: 

{zo 0.316268 + 0.129975i, .Si ^ -0.288941 + 0.157085q .Sa ^ -0.181173 - 0.05629H}. 
The corresponding periodic trajectory is shown Fig. 3 



Figure 3: Periodic Cycle solution of period 3 for a = z and /5 = 2 + 3i. 

For the parameter a = | + z and /5 = 2 + z a periodic cycle of period 10 has been 
encountered. The periodic trajectory is 

Zo -0.0197446 - 1.28723z, .si ^ 1.03398 + 0.925847z, ^ -0.406487 + 0.128166z, 

^3 ^ -0.125003-0.00142325Z ,.24 ^ 0.63328-0.516259z, ^ 0.83925+0.756558z,^ 

-0.223017 + 0.36021Z, Z 7 ^ -0.116754 + 0.0893373z, zg ^ -0.239031 + 0.16474z, zg ^ 
—0.382611 — 0.236912z. The plot of the trajectory is given in Fig.4. 



Figure 4: Periodic Cycle solution of period 3 for a = z and /? = 2 + 3z. 

Keeping fixed /? as 1 computationally it is also seen that for a = (15, 26) and (55,95), al¬ 
most all solutions corresponding to the initial values Zq, z_i G D converges to the periodic 
point (-356.366,-194.0009) and (11.6656, -0.1928) of period 55 and 199 respectively. 
Certainly there are many more such examples of a for which the same happens. The 
plot of the periodic trajectory are given the Fig. 5 and Fig. 6 for the above two examples. 
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Blue: Real Part, Green: Imaginary part 



-500 


0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

Number of Iterations x lo'' 



Figure 5: Periodic Cycle solution for a = (15,26) and 13 = 1. 


Blue: Real Part, Green; Imaginary part 



Number of Iterations x 10“ 


10000 



Figure 6: Periodic Cycle solution for a = (55,95) and /? = 1. 


In both the Fig. 5 and Fig. 6, the two figures are shown with different number of iterations 
of the periodic trajectory. 

Now we shall demonstrate few computational examples where periodic cycle exists with 
very high period roughly of order 1000. We choose a = (35, 94) and (3 = (88, 55) and for 
arbitrary complex initial values Zq and z_i, the trajectory of periodic cycles are plotted 
and corresponding phase space also plotted in the Fig. 7. 


10 





































































































100 



150 

100 



-100 ■ 

-150 I-'-'-'-'-■- 

-150 -100 -50 0 50 100 150 

Imaginary 






Imaginary 

Figure 7: Periodic Cycle solution for a = (35, 94) and (5 = (88, 55). In Each figure 5 set 
of initial values are taken and corresponding trajectories are plotted. 
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From the computational aspect, in each of the plots it is evident that for a = (35, 94) 
and P = (88, 55) the delay logistic difference equation possess very high order periodic 
cycles eventually for almost all initial values. 

4 Chaotic Solutions 

Finding chaotic solutions for the delay logistic equation 0 is interesting indeed since 
in case of real parameter a and P and the initial values Zq and there does not exists 
any chaotic solutions. 

The method of Lyapunov characteristic exponents serves as a useful tool to quantify 
chaos. Specifically Lyapunov exponents measure the rates of convergence or divergence 
of nearby trajectories. Negative Lyapunov exponents indicate convergence, while pos¬ 
itive Lyapunov exponents demonstrate divergence and chaos. The magnitude of the 
Lyapunov exponent is an indicator of the time scale on which chaotic behavior can be 
predicted or transients decay for the positive and negative exponent cases respectively. 
In this present study, the largest Lyapunov exponent is calculated for a given solution 
of finite length numerically m- 

We are looking for complex parameter a and P for which for every initial values the 
solutions are chaotic. If we consider the parameter P = 1 then the difference equation 
is known as Pielou’s equation and which is well studied earlier in case of real line. Here 
are few examples which we came across computationally. 


Parameter a, /3 = 1 

Internal of Lyapunav exponent 

a = (8,43), p = 1 

(1.205,2.623) 

a = (1,97), P = 1 

(1.845,3.028) 

a = (6,53), P = 1 

(0.785,1.718) 

a = (12,50), P = 1 

(0.373,1.485) 


Table 1: Cycle solutions of the Pielou’s equation (/3 = 1) for different choice of a and 
initial values. 


The largest Lyapunav exponents of the solutions for different initial values are lying in 
the positive intervals as stated above in the table. This ensures that the solutions are 
chaotic. It is observed that the chaos is bounded. 
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Number of Iterations x 10“ 


Figure 8: Chaotic Solutions for the Pielou’s equation (/3 = 1) of four different cases as 
stated in table 1. 
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5 Some Interesting Nontrivial Problems 

The computational experiment endorses us to pose the following important open prob¬ 
lems in this context. 

Open Problem 5.1. Find out the subset of the D of all possible initial values zq and 
Zi for which the solutions of the delay logistic equation possess chaotic solutions for a 
given parameter a and (3. 

Open Problem 5 . 2 . Find out the complex parameters a and ft such that for any initial 
values Zq and Z-i from the D the solutions of the delay logistic equation are chaotic. 

Open Problem 5.3. Find out the parameters a and jd such that for any initial values 
Zq and Z-i from the D the solution of the difference equation are periodic (globally). 

Open Problem 5.4. Characterize the parameters a and ft such that for any initial 
values Zq and z_i from the D the solution of the difference equation are periodic (globally). 
How large the period could be? Is it possible to find an upper bound? 


6 Future Endeavours 


In continuation of the present work for a generalization of the delay logistic equation, 
with varies a and /5, where I and k are delay terms and it demands similar 

l“rp.2^n —fe 

analysis which we plan to pursue in near future. 
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